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Abstract — We have made high-resolution x-ray micro- 
calorimcters using superconducting MoAu bilayers and Nb 
meander coils. The temperature sensor is a Magnetic Penetration 
Thermometer (MPT). Operation is similar to metallic magnetic 
calorimeters, but instead of the magnetic susceptibility of a 
paramagnetic alloy, we use the diamagnetic response of the 
superconducting MoAu to sense temperature changes in an x-ray 
absorber. Flux-temperature responsivity can be large for small 
sensor heat capacity, with enough dynamic range for 
applications. We find models of observed flux-temperature 
curves require several effects to explain flux penetration or 
expulsion in the microscopic devices. The superconductor is non- 
local, with large coherence length and weak pinning of flux. At 
lowest temperatures, behavior is dominated by screening 
currents that vary as a result of the temperature dependence of 
the magnetic penetration depth, modified by the effect of the non- 
uniformity of the applied field occurring on a scale comparable to 
the coherence length. In the temperature regime where 
responsivity is greatest, spatial variations in the order parameter 
become important: both local variations as flux enters/leaves the 
film and an intermediate state is formed, and globally as 
changing stability of the electrical circuit creates a Meissner 
transition and flux is expelled/penetrates to minimize free energy. 

Index Terms — About four key words or phrases in 
alphabetical order, separated by commas. 


I. INTRODUCTION 

I T HAS long been realized that a sensitive thermometer can be 
built based upon the temperature dependence of the 
magnetic penetration into a superconductor. Many groups 
have explored this concept [1-5], but difficulties in obtaining 
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desirable properties in thin film superconductors have 
prevented significant experimental success. Instead, groups 
interested in this approach have tended towards developing 
what are known as microwave kinetic inductance detectors 
(MKIDs), that are similar, but are inherently less sensitive to 
equilibrium temperature changes (this is due to their operation 
at temperatures far below the superconducting transition 
temperature and use of very high frequency read-out). Instead 
they utilize a non-equilibrium scheme with good sensitivity, 
and excellent multiplexing potential, but which arguably [6] is 
not conducive to achieving the very highest levels of spectral 
resolution. 

One of the biggest difficulties in designing a high- 
sensitivity Magnetic Penetration Thermometer (MPT) has 
been controlling the width of the temperature range of the 
superconducting transition. We have reported a solution to this 
long-standing challenge that gives a transition width suitable 
for producing a detector with usable dynamic range [7-8]. We 
use a MoAu bilayer, with a Tc of 140 mK, as our 
superconducting temperature sensor, just as is used in many 
transition-edge-sensor microcalorimeters. By subjecting the 
superconducting sensor film to the non-uniform magnetic field 
from the meander, we broaden the superconducting-to-normal 
transition enough for our device to perform as a broadband x- 
ray detector. Interestingly this broadening has not been as 
large as might be expected from the inhomogeneous field, due 
to the long-range coherence within the superconductor, so that 
the temperature dependence remains strong enough for use as 
an x-ray microcalorimeter in which 6 keV x-ray events 
typically produce temperature rises of ^ 2 mK. 

We have described our measurements and analysis of the 
transition width in our MPT microcalorimeters in a conference 
presentation [9]. Here we present details of that analysis, 
along with new results using the electromagnetic response 
kernel of BCS theory to model the inductance of the devices, 
in which the applied magnetic field varies on a scale similar to 
the superconducting coherence length. 

n. MPT Test Device and Measurement Apparatus 

The design and operation of our MPT microcalorimeters are 
similar to Metallic Magnetic Calorimeters (MMCs) [10]. In 
both MPTs and MMCs we use the temperature dependence of 
the magnetic susceptibility of a metal film to change the 
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inductance of a superconducting niobium meander coil 
coupled to a SQUID sensor. However, instead of the 
paramagnetic susceptibility of a Au:Er alloy, we use the 
diamagnetic susceptibility of a superconducting MoAu bilayer 
to sense temperature change in an overhanging x-ray absorber 
in electrical contact with the sensor film. 

Figure 1 shows one of the MPT devices described here. We 
used a planar niobium meander with 2.5 pm line width, 5.0 
pm pitch, and 0.45 pm thickness. We covered the meander 
with a 0.32 pm oxide insulation layer, and then a MoAu 
bilayer film 0.38 pm thick. The bilayer had a superconducting 
transition temperature T c = 140 mK, and a residual resistance 
ratio RRR =8.1, with its normal state resistivity dominated by 
the high conductivity of its 0.33 pm thick Au component. 
Details of the fabrication methods are similar to what we 
employ for making MMCs [11]. 



Fig. 1. Magnetic Penetration Thermometer without x-ray absorber. A thin- 
film of MoAu bilayer covers a niobium meander coil, with 5 nm pitch, in a 
100 pm x 100 pm region. An aluminum oxide layer insulates the coil from the 
bilayer. Niobium microstrip lines from SQUID input and MPT bias circuit 
(top and right) supply a current to the meander coil with a via through the 
oxide. Th: current from the meander returns (at lower left) to the ground of 
the micros trip in a split path, symmetrically along the top and bottom edges. 

We mounted in the MPT devices in a magnetically shielded 
package in an adiabatic demagnetization refrigerator. With 
short wire bonds designed to minimize stray inductance, we 
connected the MPTs to two-stage dc-SQUID amplifiers from 
PTB-Berlin [12]. At a temperature > 1 K, we stored a 
persistent current, ranging up to several milliamps, in the 
niobium circuit biasing a series array of MPTs, then cooled to 
temperatures as low as 0.030 K. We measured the dc level and 
noise spectra of the flux <p in the SQUID input at various bias 
and temperature combinations, and recorded pulses in 
response to applied x-rays with energies of 1.5 or 5.9 keV. 
Further details of the measurement apparatus are given in [7]. 

m. Experimental Results 

During initial cooling after storing a current in the niobium 
bias circuit, we observed a staircase of flux jumps in the 
SQUID dc level that was repeatable if the temperature was 


cycled above 140 mK. However, for bias currents < 1.25 mA, 
the flux-temperature curve was smooth upon warming, and 
largely non-hysteretic thereafter for temperature excursions 
below 140 mK. At a stored bias current 7 b = 1.25 mA, the 
maximal flux change was 14 <P Q . The greatest responsivity to 
temperature was « 1000 <PJ K, which occurred near T = 80 
mK. At higher bias currents the peak flux signals decreased. 
For three different values of the bias current stored in the 
niobium circuit, the bold curves in Fig. 2 show the SQUID 
input flux change from high temperature, normalized by bias 
current, as the MPT was warmed. 
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Fig. 2. Measured flux, per unit of bias current, plotted against temperature 
(bold curves) for stored bias current values 0.50, 1.00, and 1.25 mA. Also 
shown are various model curves, discussed in Sec. IV. 

With the bias circuit wired in a different configuration, we 
were also able, during a separate refrigerator run, to apply the 
bias current first at 30 mK and then measure the SQUID flux 
change during warming to T > 140 mK, instead of our normal 
procedure of cooling after storing the bias current with the 
MoAu film in the fully normal state at T > 1 K. Comparison of 
the two total flux change values indicated that > 80% of the 
flux penetrating the normal bilayer is expelled for 7 b = 1.25 
mA during our normal bias procedure. 

Noise measurements [7] with at 7 b = 1.25 mA and T = 80 
mK showed the presence of no significant excess noise, only 
SQUID and phonon noise. The phonon noise component 
extended up to 5 kHz, and was consistent with the expected 
electron-phonon and Kapitza thermal conductance between 
the bilayer and substrate. The measured ratio of heat capacity 
to temperature ~1.5 pJ/K 2 was similar to the normal state heat 
capacity of the gold contained in the bilayer, but measurement 
precision was limited by the accuracy of our measurement of 
d<p/dT at the dc operating point during the noise measurement. 

Noise spectra at 4 K with the niobium circuit 
superconducting but the aluminum wire bonds normal gave 
the total circuit inductance at that temperature to be 4.14 nH. 

IV. Modeling Methods and Results 

A. Circuit Model 

From quantization of the flux in the niobium meander 
coil/SQUID input circuit loop, and the near constant bias 
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current supplied to the parallel combination of the SQUID 
input coil and the MPT meander (enforced by a large ballast 
inducto; in the bias circuit), one can derive the following 
relations between the meander coil inductance L m (T) and the 
currents /* and I m in the SQUTD input and meander: 


/. = 


-UA>)-Ur)] 


L m (T) + L, 

! =7 7 

m h " L m (T) + Li 


0 ) 

( 2 ) 


The SQUID input coupling was Mi = l # 0 per 5.4 pA, and 
its input inductance L* was known to be 1.8 nH. The 
inductance of the meander with the MoAu fully normal was 
computed to be 0.76 nH. From 4.14 nH total circuit inductance 
at 4 K. the value of stray inductance from wirebonds and 
microstrip was therefore 1.58 nH. The total input inductance 
including strays was L v = 3.38 nH. The overall change in <p = - 
MJ i between hot and cold limits in Fig. 1 constrains L m (0) = 
0.55 nH. Figure 3 shows the experimental meander inductance 
versus temperature in comparison with models. 


B. Screening Currents 

For 1 j < 1.25 mA, we observe a regime at the lowest 
temperatures, T < 60 mK, where the MoAu seems to not trap 
the field applied by the meander, but to expel the flux as 
completely as possible by means of screening currents. In that 
regime, the SQUID flux is linear in the stored bias current, i.e. 
<j>/ / b is nearly independent of 7 b . 

Measurements of the properties of Transition Edge Sensors 
made from similar MoAu bilayers have shown that such films 
exhibit non-local superconducting behavior, with a coherence 
length §) ~ 0.7 pm much larger than their London penetration 
depth ^ ~ 0.08 pm [13]. If we instead treat our MPT bilayer 
as a local superconductor, we can easily compute the expected 
meander inductance resulting from screening currents in the 
bilayer by solving the London equation numerically using the 
method of Chang [14] (see Fig. 3a). However, we then predict 
much less variation of the SQUID flux and meander 
inductance at the low-T end of our data than we observe in 
Fig. 1 and Fig. 3b. 
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Fig. 3. (a) Example of computed current densities on Nb and MoAu film 
cross-sectbns for London model, (b) Measured meander inductance versus T 
and fitted curves for London (local limit) and BCS (non-local) models. 


We have computed the expected meander inductance in the 
non-local limit using the electromagnetic response kemal from 
BCS theory. We consider a superconducting slab of finite 
extent in the z direction and infinite extent in the x-y plane. 


There is a sheet current pointing in the y direction with spatial 
variation in the x direction offset from the superconducting 
slab in the negative z direction. The current sheet is 
representative of a meander coil underneath a superconducting 
body. Our goal is to compute the inductance associated with 
the current sheet in the presence of the superconducting slab. 
To do this, we will compute the magnetic field in the region 
and the supercurrcnt in the slab and equate the total energy 
with the energy stored in an inductor. 

We start with two equations, Ampere’s law neglecting the 
displacement current and the BCS expression for the 
supercurrent as the convolution of the BCS integration kernel 
with the electromagnetic vector potential. The current in 
Ampere’s law is the sum of the sheet current and the 
supercurrent. Due to the form of the supercurrent, it is 
mathematically most convenient to work in terms of the vector 
potential. This is done by replacing the magnetic field in 
Ampere’s law with the curl of the vector potential. We choose 
to work in the Coulomb gauge so that the curl of the magnetic 
field becomes the Laplacian acting on the vector potential. 
The independence of the current on the y direction allows for a 
simplification of Ampere’s law. We can satisfy the gauge 
condition by allowing only the y component of the gauge field 
to be non zero and such that it is a function of only the x and z 
variables. A further simplification can be arranged by 
expanding the current and vector potential in a Fourier series 
in x coordinate. This reduces the differential equation for the 
vector potential to a one-dimensional equation for each of the 
z-dependent Fourier modes, with the Laplacian becoming a 
one-dimensional modified Helmholtz operator. The 
Helmholtz operator is inverted using its Green’s function 
converting the differential equation into an integral equation. 
The integral equation is discretized and the solution for the 
vector potential is found numerically. The solution is used to 
construct the magnetic field profile and the supercurrent 
profile, which is used to compute the inductance of the setup 
as stated above. 

The model has three parameters, the superconducting gap of 
the slab, the Fermi momentum, and quasiparticle effective 
mass of the superconducting material. The Fermi momentum 
and quasiparticle effective mass can be expressed in term of 
the zero temperature London penetration depth and the 
coherence length, which are more convenient parameters to 
use. In the absence of independent means for determining the 
zero temperature London penetration depth and coherence 
length, the two parameters are fit to the data. 

C. Meissner Effect and Landau Model with Bias Circuit 

At temperatures above 60 mK, we observe a bias-current- 
dependent transition temperature at which the reduction of the 
meander inductance by effects in the MoAu bilayer rapidly 
decreases. We interpret this as a manifestation of the Meissner 
effect specific to the geometry and bias circuit for our devices. 

F. London analyzed [15] the Meissner effect for a thin film 
placed in a uniform applied field in terms of the free energies 
of the fully normal state and the superconducting state with 
screening currents (computed in the local limit). Applying 
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London’s method to our films would predict sudden 
transitic ns (above a critical temperature varying with / b ) to <p = 
0 from a universal low-T curve for (p / 7 b . Our observed 
transitions are softened on both the warm and cold sides. 

The transition would be modified if the MoAu energy gap is 
reduced in response to the magnetic field. As investigated by 
Douglass [16-17], a simplified Ginzburg-Landau model of this 
effect predicts that the character of the phase transition 
changes to second-order if the field increases to the critical 
value with a sufficiently small ratio of film thickness to 
effective penetration depth. Curves based on the G-L model 
used by Douglass, ignoring the spatial variation of field, are 
compared with our data in Fig. 1. The reduction of the energy 
gap by the field does soften the flux-temperature curves on the 
low-T side, but not as much as observed, while on the high-T 
side the predicted transition is infinitely sharp. 

We obtained much better qualitative agreement with the 
observed transitions using a Landau phase transition model. 
We assumed that the Type-I MoAu bilayer can enter an 
intermediate state in which stripes of bilayer go normal (gap 
drops to zero over a coherence length) as shown in Fig. 4a. 
When coupled to the bias circuit, we find the free energy 
difference between the intermediate state and the fully normal 
state is given by, 






w'here the terms, in order are: the free energy if all the MoAu 
were fully normal, the magnetic field energy of the meander if 
the MoAu were fully normal, the extra magnetic and kinetic 
energy for a fraction (\-j) of superconducting MoAu screening 
the Nb coil, the negative of the condensation energy of the 
superconducting volume of MoAu, and the G-L interface 
energy between normal and superconducting regions. The 
length of the sides of the MoAu square are s , while d is the 
MoAu thickness and Z? c (7) is the thermodynamic critical field. 
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Fig. 4. (a) Postulated intermediate state in MoAu is high-lighted in sketch of 
Nb meander coil, (b) Free energy curves at 1.25 mA bias, for temperatures 
increasing from 77 to 8 1 .9 mK in 0.7 mK increments, plotted versus fraction / 
of coil in normal state. Location of minimum free energy is shown by a cross, 
and range of f values expected from thermal fluctuations is shown by circles. 

Our Landau model of the Meissner transition predicts, with 
reasonable parameters, a finite maximum value for d^dTthat 
is of a similar magnitude to the observations at high bias 
currents (and T substantially less than the zero-field Tc = 140 
mK). The model requires that the edge barrier for fluxon 


escape from the MoAu must be small and that pinning of the 
flux in the bilayer is weak. 

V. Conclusion 

MPTs enable x-ray microcalorimcters with zero bias power 
dissipation and excellent energy resolution. MPTs can be 
optimized to provide dtp/dT as large as 1000 tfVK, with no 
excess noise, thereby reducing the importance of SQUID 
noise, provided the sensor film is a Type-I superconducting 
film with weak flux pinning. The temperature region of 
steepest dtpIdT represents the Meissner effect in the small 
device as flux is expelled/penetrates to minimize free energy. 
The transition width, and microcalorimeter dynamic range, 
can be increased by non-local superconducting effects when 
the coherence length is similar to the scale of spatial variations 
in the field applied by the meander coil. 
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